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We investigate the quantum dynamics of a periodically kicked Bose-Einstein Condensate confined 
in a one dimensional (ID) Box both numerically and theoretically, emphasizing on the phenomena 
of quantum resonance and anti-resonance. The quantum resonant behavior of BEC is different from 
the single particle case but the anti-resonance condition (T = 2n and a = 0) is not affected by 
the atomic interaction. For the anti-resonance case, the nonlinearity (atom interaction) causes the 
transition between oscillation and quantum beating. For the quantum resonance case, because of 
the coherence of BEC, the energy increase is oscillating and the rate is dramatically affected by the 
many-body interaction. We also discuss the relation between the quantum resonant behavior and 
the KAM or non-KAM property of the corresponding classical system. 
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I. INTRODUCTION 

Quantum systems under a periodically driving force 
are of great interest in varied fields of physics for its 
versatile applications in microscopic manipulations and 
control lj. Its dynamics demonstrate many interesting 
behaviors, like dynamical localization and chaos-assisted 
tunneling, to name a few 0. [J [IH3 Among them, quan- 
tum resonance (QR) and anti-resonance (AR) are two in- 
teresting phenomena 0,0- QR says that under a certain 
resonance condition, a particle acquires energy from an 
external force most efficiently leading to its energy in- 
crease with time in a square law. In the other limiting 
case, the AR case, the particle will bounce between two 
states and its energy shows a periodic oscillation. 

QR and AR are pure quantum behaviors without clas- 
sical counterparts. In the well-known kicked rotor sys- 
tem, given a value of the kick strength K, special res- 
onant regimes of motion appear for periods with values 
T = 47t| , wnere t ne integers p and q are mutually prime. 
Under these conditions the system regularly accumulates 
energy which grows quadratically in the time asymptotic 
0- The case f = 3 presents a completely periodic be- 
havior with period 2T. This is the AR case. QR and 
AR has been observed in the atom optics imitation of 
the quantum kicked rotor 0. Recent realizations of the 
BEC jfj make us curious about whether or not QR and 
AR also exist in BEC systems and how the non-linearity, 
stemming from mean field treatment of the atomic inter- 
action 0| , affects quantum resonances. Recently in 
it has been shown that a BEC in a quasi ID box can be 
achieved. This experiment provides a good condition to 
investigate the quantum resonances of BEC. The S kick 
can be realized using counter propagating laser beams 



and its spatial shape can be adjusted by phase mismatch 
of laser beams. The interaction strength between atoms 
g can be changed using a Feshbach resonance. These mo- 
tivate us to study the quantum resonances of BEC under 
this experimental condition. 

In this paper, we consider a BEC trapped in a ID 
box and kicked periodically, and study how the atomic 
interaction affects the quantum resonant behaviors. We 
find that the QR and AR conditions for this system are 
different from the quantum kicked rotor system. For this 
system the AR can be only found for a special spatial 
shape of the S kick when T = 2ir (i.e. 2 = i), for other 
shapes of kick and T = 47r| 2n) the QR will be 
observed. Because of the coherence of BEC, for the 
AR case, the nonlinearity (atom interaction) causes the 
transition between oscillation and quantum beating; for 
the QR case the energy increase is oscillating and the rate 
is dramatically affected by the many-body interaction. 
We also find that for the QR case the nonlinearity (atom 
interaction) suppress the sensitivity to the spatial shape 
of the kick. Finally, we discuss the relation between the 
quantum resonant behavior and the KAM or non-KAM 
property of the corresponding classical system. 

This paper is organized as follows. In section II, we 
introduce the model. In section III, we show how the in- 
teraction between atoms in a BEC changes the evolution 
of the energy in quantum resonant cases. In section IV we 
present our analytical results which show how the many- 
body interaction affects the evolution of the energy. In 
this section, we also show why for certain values of the 
parameter a there is no AR and we show that this is not 
related to whether the corresponding classical system is 
KAM or non-KAM. Finally, in section V, we present our 
conclusions. 



2 



II. THE MODEL 

We focus our attention on the dynamics of a quasi- ID 
BEC confined in a cigar shaped trap with a pulsating 
potential. In the limit of validity of the mean-field treat- 
ment, this system can be described by the dimensionless 
non-linear GPE: 



H = - 



2dx 2 



+ g\<j>(x,t)\ 2 + Kco$(x + a)5 T + U(x) (1) 



where U(x) = oo, for x < 0, x > n, and U(x) — else- 
where, g — aiD^h Na/m is the scaled strength of non- 
linear interaction, N is the number of atoms, a is the 
s-wave scattering length, a±o is a coefficient which com- 
pensates for the loss of two dimensions K is the 
kick strength, 8 t (T) represents ^2 n 5(t — nT), T is the 
kick period, and x denotes the position on the x axis. 
The variable x £ [0, tt] and a is a parameter between 
and 2n which can be controlled by the phase mismatch 
of laser beams. Due to symmetries, the only important 
interval for the parameter a is a £ [— tt/2, tt/2]. Because 
the BEC is unstable under kicks if the nonlinear interac- 
tion is large [l3| , we only study g with a maximum value 
of 0.5. This value is very likely in the stable region where 
the number of condensed particles is much bigger than 
the number of non-condensed ones for a long enough time 
that we can use the GPE to study the evolution of the 
wave- function. 

In the resonant case the energy grows in average 
quadratically and at the same time the number of non- 
condensed particle is also growing. This means that our 
results are valid for a limited number of kicks because af- 
ter that we lose the coherence of the condensate. It could 
be interesting to see experimentally where this limit is. 
The energy of a particle is given by (E) — 

J* dx[(j>*(— ^2 + ||0| 2 )<^]. The evolution of the wave- 
function is given by numerical integration of Eq.Q, 
over a certain number of kicks using the split-operator 
method. 

In our study, we use the ground state of the Hamilto- 
nian H = — ^fj? + 9\4>( x > t)\ 2 (with the same boundary 
condition given above) as the initial condition. Due to 
the shape of the potential U(x), <fi(0,t) — 4>(~K,t) = 0. 
The wave function <fi(x, t) satisfies the normalization con- 
dition J* \cj)\ 2 dx = 1. For a positive g the ground state 
is given by (01): 



m (2K(m) 
9 V n 



(2K(m)-\i 



(2) 



where K(m) is the elliptic complete integral of the first 
kind and sn(x, m) is the Jacobi elliptic function. The 
parameter m is included in the interval [0, l]and is related 
to g by: 



— (2K(m)) 2 
gn 



1 - 



E(m) 
K{m) 



= 1. 



(3) 



which comes from normalization condition. For negative 
values of g, the initial condition is given bv|l4|: 

m f2K(m)\ ( , w 2x \ 
- (— [ -2Jcn(K(m)(— -l)\m) (4) 

where cn(x, m) is the elliptic Jacobi function an and m 
and G, which is G = —g, are related by: 



(1 - m) = 1. (5) 



III. QUANTUM RESONANCE AND 
ANTI-RESONANCE OF BEC 

A. Anti resonance 

For this model we find that the AR can only be ob- 
served with the condition T = 2ix and a = (shown in 
FigUJ. This condition is different from the one of kicked 
rotor studied by Zhang et al where they discovered 
that the AR condition is for T — 2tt and independent 
on the shape of kicks. In FigQ] one can see that if the 
non-linear term is zero, the energy oscillates in time with 
a period 2T. However, when non-linear term is nonzero, 
the energy oscillates in time with an amplitude that de- 
creases gradually to zero and then revives, similar to the 
phenomena of beating in classical waves. 

For the quantum beating case there are two frequen- 
cies: one is the frequency of kick and another is the beat- 
ing frequency which is due to the coherence of the BEC 
and can be obtained approximately by a 2-state model. 

As for the one in 13], for a = 0, our model can be 
mapped onto a 2-state model (0], 0], H3)- We can 
write the wavefunction as a sum of only the ground state 
and the first excited state with relative population a and 
b (with normalization condition \a\ 2 + \b\ 2 — 1). Defin- 
ing S z = \a\ 2 — \b\ 2 as the population difference between 
ground state and first excited state and — arctan^y/S'z) 
as the relative phase between the two states, we can ex- 
press the Hamiltonian as: 



H-- 



:S Z 



9_ 

2n 



4 



K 



S t (T)S x 



(6) 



Then, one can obtain the beating frequency of the evo- 
lution of the energy from the 2-state model, 



fbeat 



gcos(K) 



(7) 



We can see in Fig. |21 that the theoretical approximate 
agrees very well with numerical calculations. 

The quantum resonant behaviors of kicked BEC in ID 
box are different from the one of kicked rotor studied 
by Zhang et al 13]. For the kicked BEC in ID box the 
QR behaviors can also be controlled by the spatial shape 
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FIG. 2: (a) Modulation frequency versus kick intensity (K) for 
fixed interaction value (g = 0.2). (b) Modulation frequency 
versus g for fixed kick-intensity (K — 0.5). The theoretical 
result of Eq.0 agrees very well with the numerical simula- 
tion. 



TIME (in unit of kicks) 

FIG. 1: Energy ((E)) evolution for different values of the 
interaction strength g with fixed K = 0.5. The period is fixed 
to T = 2n and a — which corresponds to AR condition. 
For g = (a) the evolution of energy is perfectly periodic 
with period 2. For the other two cases (g — 0.1 in (b) and 
g — 0.3 in (c)), the evolution is quasi-periodic and we can see 
the phenomenon of beating. 

of the kick which can be adjusted by mismatch of laser 
beams. In our model the spatial shape is parameterized 
by a. Following we will show that if a ^ 0, we will have 
QR even if T = 2tt (AR condition of kicked rotor). 



B. Quantum resonance 

For any nonzero a, the quantum behavior of the system 
is very different from the case of a — 0. In Fig. [3] we 
show the energy evolution with time (in unit of kicks) 
for different g = (a), g = 0.1 (b) and g = 0.3 (c) with 
a = 0.1. We can see that for g — 0, AR does not exist 
anymore and there is only QR. The energy increase with 
time on average in square law. 

For j / 0, because of the coherence of BEC, the 
energy increase is oscillating and the rate is dramatically 
affected by the interaction term. However, though the 
energy oscillates, the energy on average has a quadratic 
increase as for the resonant case. 

The behavior of the energy with the number of kicks 
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FIG. 3: Energy ((E)) evolution for K = 0.5 and a = 0.1. 
The interaction strength is g = for (a), g = 0.1 for (b) and 
g — 0.3 for (c). The motion is neither periodic nor quasi- 
periodic. In (b) and (c) we can also see the oscillation due to 
the non-linear interaction. 



for different values of a and g is summarized in Fig. 0] 
where the value of energy after 50 kicks for K — 0.5 and 
different values of g and a is shown. We can see many 
interesting things in this figure. First of all, there is a 
symmetry axis at a = n. For g = the value of energy is 
symmetric and the symmetry is broken by the non-linear 
term. Opposite values of g are symmetric (considering 
only the kinetic energy) with respect to this axis. The 
energy reaches a maximum at a — tt/2 for negative g 



and at a = 3n/2 for positive g. The breaking of the 
symmetry is due to the coherence of the BEC. 
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FIG. 4: Energy ((E)) evolution for different values of g with 
fixed K — 0.5. The interaction strength is g = (blue), 
g = ±0.1 (light blue and green), g = ±0.3 (red and purple) 
and g = ±0.5 (black and yellow). We can see a symmetric 
behavior for opposite values of g as long as the strength is 
not too large (case \g\ = 0.5). Positive values of g reach a 
maximum for a — 3n/2 and negative for a — tt/2. 

From Fig. QJit seems that, the bigger the value of g, the 
lesser is the effect of a small change of a. We can see this 
in detail in Fig. [5] where we show for g — 0, g = 0.1 and 
g = 0.3 the effect of small values of a. For a — 0.4 only 
the case with g = 0.3 is still stable, while for example for 
g — 0, a tiny perturbation such as a = 0.01 is enough to 
make the energy increase rapidly. This shows that, the 
bigger is the interaction between the atoms in the BEC, 
the less sensitive to a variation of a is the system. 



IV. ANALYTICAL STUDY 

A. Approximate study of the influence of the 
interaction 

In this section, we shall demonstrate that it is possi- 
ble to understand analytically why the evolution of the 
energy has the above behavior. As it will be shown, the 
underlying mechanism is the interaction of the atoms in 
the BEC. 

We approach the problem by looking for an expression 
of the energy to the first order in g. For small values of g 
we can approximate the shape of the wave-function (f>(x) 
by: 
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FIG. 5: Energy ((E)) as function of the number of kicks for 
different values of g and a = (black), a = 0.02 (red) and 
a — 0.04 (green). The interaction strength is g — (c), 
g = 0.1 (b) and g = 0.3 (a). The bigger is g the more stable 
is the motion to changes of the parameter a. 



Mapping this model on a periodic ring, we can predict 
the value of the energy after a certain number of kicks. 
In general a model with x £ [0, w] and an infinite square 
potential can be mapped into a model with x £ [0, 2-7r] 
and periodic boundary conditions. The new initial 
condition (f>(x,0) is, for < x < n, 4>(x,0) — (j)(x,0)/y/2, 
and for n < x < 2n, <f>(x,0) = <j>{2ir - x,0)/s/2, 
where <j){x, 0) is the initial condition for the model with 
x € [0, 7r]. The new kick V(x) is given, for < x < 7T, by 
V(x) = V(x) and, for tt < x < 2tt by V(x) = V(2ir - x) 
where V(x) in this case is K cos(x + a). This mapping 
could be done for any potential. 

The evolution of the new wave-function, after an even 
number (N = 2M) kicks, is given by: 



(x, NT) = e 



-iM(V+V„) 



1 ± e~ i2 f sin2 ( x ) 



N 



4>(x,o) 



(x,T) 



= e -^) e -*9H 2 i e ^T e -i ff |0| 2 2: 



(9) 

■ 0) (8) where V v = V(x + n). In this way we can compute the 
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energy which is given by: 

, . N 2 K 2 „ s 16 1 3 , N 

E(N) = — — sin 2 (a) - —gsma + - + —g 10 

8 157T 2 107T 



To obtain this result, we have used <f)(x, 0) = 
which is a very good approximation in the case of small 
values of g here studied. 

Eq. HI Op approximates the value of the energy after an 
even number, N, of kicks. We can see in Fig. |S|how good 
the approximation is. From IjlOp we can see if a = the 
energy is independent of the number of kicks N and if 
a / the energy increases with the number of kicks in 
square law. 

Moreover it is possible to see in Eq. (|10p why for g = 
the behavior is completely symmetric and also where the 
symmetry for opposite values of g comes from. However, 
it is not possible to understand the behavior for bigger 
values of g in this way. For example we can see in Fig. 
0] that for g — ±0.5 the maximum is smaller than for 
g = ±0.3. It is also not possible to show the oscillation 
of the energy due to the non-linear term, but for our 
purpose of understanding the average growth of g for 
different values of a, the method shown in this paper is 
good, at least for small values of g. 



B. Classical dynamical properties and quantum 
resonance 

The study of quantum systems whose classical coun- 
terpart is non-KAM already showed interesting results. 
In both the classical and quantum cases, the non-KAM 
systems 0, 0, l^jj, demonstrate quite different be- 
havior from the KAM one. For instance, in classical 
KAM systems, as the external or driven parameter is 
increased, the invariant curves gradually break up. Lo- 
cal chaos evolves into global chaos and diffusion takes 
place. In a non-KAM system, there are no such invari- 
ant curves for any small external or driven parameters. 
Quantum mechanically, the quantum interference sup- 
presses the classical diffusion leading to a so-called expo- 
nential localization^, while in a non-KAM system, the 
localization becomes power law, or in other words, there 

is no localization^]- 

In our case, it is interesting to notice that the system 
that we study presents AR when the classical counterpart 
is KAM (for a = 0), and that there is QR when the 
calssical equivalent system is non-KAM (for a ^ 0). We 
would like to understand whether the properties of the 
quantum system and its classical equivalent (AR with 
KAM and QR with non-KAM) are related. 

To understand this we start showing a general result 
which is valid for a generic periodically kicked systems 
with x £ [0, 2tt] and the kick is given by a generic V(x). 
Now, starting from Q it is easy to see that 4>(x)n+i = 
cxp(— iV(x))<f)(x + 7t)at. So after two kicks with a period 
T = 2ir we have that: 
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FIG. 6: Evolution of the energy ((E)) for different values of 
the interaction strength g and K — 0.5 after 50 kicks for dif- 
ferent values of a. The numerical (black) curve is compared 
to the analytical ones given by Eq.JTTj| (red). The red curve 
approximates the numerical results. The values of the inter- 
action strength are g = 0.01 (a), g = —0.01 (b), g = 0.05 (c), 
g =-0.05 (d), = 0.1 (e), <? = -0.1 (f). 



<j)(x) N+2 = exp(-i(V(x) + V(x + 7t)))<p(x) N (11) 

From this point it is easy to derive that the condition 
for AR is T = 2tt and 



V(x) + V(x + tt) = C 



(12) 



This result, after doing the mapping discussed in Section 
III, can be applied to our model which, as it should 
be emphasized, is characterized by an infinite well 
and x £ [0, 7r]. For a ^ or a ^ ir, the mapping 
on the circle of the kick, V, does not follow Eci. (|12p . 
V(x) + V(x + tt) 7^ C, but it is a function of x. 

It is now obvious that AR is not related to whether 
the corresponding classical system is KAM or non-KAM, 
but due to Ea. l|12p . We can have KAM systems without 
AR and non-KAM systems with AR. 

This is shown also in Fig. [7| where the kick is given ei- 
ther by V(x) = K cos(2a;) (KAM) and we see QR, either 
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by V(x) = tt/2 - x for < x < ir and V(x) = -3tt/2 + x 
for 7r < x < 2ir and we have AR. 



A 
HI 

V 



25-| 
20 
15 
10 

5 





(a) 



10 



20 



— I ' 1 — 

30 40 



50 



A 
HI 

V 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 




(b) 



10 20 30 40 
TIME (in unit of kicks) 



50 



FIG. 7: (a) Energy ((E)) evolution for a kick of type V(x) — 
Kcos(2x). T = An\ and K = 0.1. We do not have AR even 
though the system is KAM. (b) Evolution of the energy ((E)) 
with the number of kicks for a kick of the kind V{x) — n/2 — x 
for < x < 7T and V(x) — — 3-k/2 + x for tv < x < 2tt. The 
system is non-KAM but we have AR. 



V. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have investigated the quantum dy- 
namics of a periodically kicked Bose-Einstein Condensate 



confined in a ID box both numerically and theoretically, 
emphasizing on the phenomena of QR and AR. We find 
that the atomic interaction does not affect the AR condi- 
tion (or QR condition) . However, the resonant behaviors 
of BEC is different from the single particle case. For 
the AR case, the nonlinearity (atom interaction) causes 
the transition from oscillation to quantum beating. For 
the QR case, because of the coherence of BEC, the en- 
ergy increase is oscillating and the rate is dramatically 
affected by the interaction between atoms. And, the 
rate at which the energy increases in the system depends 
on the atom interaction. The interaction breaks the 
symmetric evolution of the energy for different values of 
a around a — n. We have also found that, the stronger 
the interaction between atoms, the more stable the sys- 
tem is to small changes of a. This means that for bigger 
interaction the AR behavior will be much more stable 
to errors in the matching of the kick-generating lasers 
and the trap. We also discussed the relation between the 
quantum resonant behavior and the KAM or non-KAM 
property of the corresponding classical system. 



We would like to emphasize the fact that the system 
that we studied can be realized with current experimen- 
tal techniques. In this way with the same experimental 
set-up it could be possible to observe the phenomenon 
which we have shown, such as quantum beating, and the 
phenomenon of the destruction of AR. The phenomenon 
of quantum beating can be used to measure the value of 
the interaction of the atoms, and the breaking of AR can 
be used to see if the laser matches with the trap. 
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